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ABSTRACT 

The dynamical foundations of a disk models are described. At the heart of 
the viscous formalism of accretion disk models are correlations in the fluctuating 
components of the disk velocity, magnetic field, and gravitational potential. 
We relate these correlations to the large scale mean flow dynamics used in 
phenomenological viscous disk models. MHD turbulence readily lends itself to 
the a formalsim, but transport by self-gravity does not. Nonlocal transport is 
an intrinsic property of turbulent self-gravitating disks, which in general cannot 
be captured by an a. model. Local energy dissipation and a-like behavior can 
be re-established if the pattern speeds associated with the amplitudes of an 
azimuthal Fourier decomposition of the turbulence are everywhere close to 
the local rotation frequency. In this situation, global wave transport must be 
absent. Shearing box simulations, which employ boundary conditions forcing 
local behavior, are probably not an adequate tool for modeling the behavior of 
self-gravitating disks. As a matter of principle, it is possible that disks which 
hover near the edge of gravitational stability may behave in accord with a local 
a model, but global simulations performed to date suggest matters are not this 
simple. 

Subject headings: hydrodynamics — instabilities — turbulence 



-3- 



1. Introduction 



For many years, the principle uncertainty and greatest impediment for the development 
of accretion disk theory was an understanding of the origin of turbulent transport. In their 
classic paper, Shakura & Sunyaev (1973) made the physically reasonable and enormously 
productive ansatz that, whatever the underlying cause for its existence, the turbulent 
stress tensor scaled with the local gas pressure P. They denoted the constant of 
proportionality as a, and the "a disk" moniker has since become synonymous with the 
standard disk model. Despite the ongoing development of increasingly sophisticated large 
scale numerical models, a disk modeling still remains the central link between theory and 
observations, the cornerstone of accretion disk phenomenology. 

In the last several years, a promising candidate has emerged as the physical basis for a 
disk models. This is the magneto- rotational ( "Balbus-Hawley" ) MHD instability (Balbus 
& Hawley 1998, and references therein). A large and ever-increasing body of numerical 
simulations leaves little doubt that this instability leads to the turbulent enhancement of 
angular momentum transport within an accretion disk. What has been lacking, however, is 
a systematic explanation of how turbulence may or may not lead to the phcnomcnological 
a disk equations which have been in use for many decades. Since this approach has been 
(and continues to be) the link between accretion disk theory and observations, a better 
understanding of the a formalism is clearly desirable. The present paper seeks to fill this 
role. 

The dynamical foundations of viscous disk theory have always been somewhat fuzzy, 
and the benefits of sharpening our understanding are numerous. For example, it is 
desirable to clarify which results of a disk phenomenology are truly fundamental, and 
which have more limited domains of applicability. The nature of time-dependent turbulent 
transport could be more fully elucidated: in what sense is it equivalent to a viscous 
stress? Most interestingly, by introducing the intermediate integral scales of turbulence 
into the investigation, an enormously richer class of physical problem emerges. Classical 
a disk theory addresses mean flow (macroscopic) dynamics, subsuming all integral scale 
structure into a viscous stress tensor. In this approach, macroscopic disk structure is 
coupled directly to the dissipative scales. One cannot begin to answer such questions as 
whether the turbulence is self-maintaining, or whether the transport is local or global; 
everything is simply prescribed. Finally, there are important questions facing disk modelers 
in nonmagnetized disks, such as protostellar disks on scales larger than a few AU. Is it 
sensible to model such regions with an a viscosity? Are disks which evolve under the 
influence of self-gravity a disks? What distinguishes turbulence well-modeled by a viscosity, 
from turbulence which is not? These and related questions form the focus of this paper. 

An overview of this paper is as follows. In §2, we present a review of classical viscous 
accretion theory. Although this material is well-known, we revisit it with a renewed 
attention upon how the viscous stress appears in the angular momentum and energy fluxes. 
This becomes a benchmark for the turbulent theories discussed in §§ 3 and 4. In §3 we show 
that MHD turbulence acts very much along the lines of classical viscous theory, in both 
steady-state as well as evolutionary disk models. In §4, it is shown the turbulent transport 
arising from self-gravity is not, in general, compatible with a viscous formalism. We discuss 
the physical basis for this behavior, and show that there are hmiting cases of restrictive 
generality, which are compatible with a disk theory. Finally §5 summarizes our findings. 
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2. Preliminaries 



2.1. Classical Viscous Disk Theory 



We begin with a brief review of classical viscous disk theory (Lynden-Bell & Pringle 
1974, Pringle 1981). In Cartesian coordinates, with x, y, z represented by dummy indices i, 
j, k, the viscous stress tensor takes the form (Landau & Lifschitz 1959): 

(Jij = r]{djVi + diVj - '^SijdkVk). (1) 

We use the standard notational convention of di denoting the partial derivative with respect 
to spatial coordinate i, and summation over repeated indices is implied unless stated 
otherwise. Vi is the ith component of the velocity vector, and t] is the dynamical viscosity. 
To the extent that the fluid behaves incompressibly, we may ignore the divergence term, 
a standard and generally well-justified procedure for the class of turbulence we wish to 
consider here. The idea of viscous disk theory is to regard the effects of turbulence as 
greatly enhancing the magnitude of rj beyond its microscopic value, and doing nothing else. 

The dynamical equations are mass conservation, 

dtp + djipvj) = 0, (2) 

and the equation of motion 

p{dtVi + VjdjVi) = -dj{P5ij - o-y) - pdi^, (3) 

{6ij is the Kronecker delta function) which may also be written 

dtipvi) + djipViVj + P6ij - aij) = -pdi^, (4) 

an explicit statement of momentum conservation. (Our notation is again standard: p is 
the mass density, $ is the gravitational potential, P is the gas pressure.) At this stage, 
we assume that $ is an imposed central disk potential; self-gravity is considered in §4. 
In viscous disk theories, momentum transport — or more usefully, angular momentum 
transport — is the task of aij. If the differential rotation rate decreases with increasing 
radius, viscosity transports angular momentum outward. 

Multiplying eq. ||^ by fj, integrating terms by parts, using mass conservation, and 
finally summing over leads to a mechanical energy equation 

dt{pv'^l2 + p^) + dj{pv'^Vj/2 + p$ + Pvj - Viaij) = PdjVj - {djVi)aij, (5) 

where 

v'^ = ViVi (6) 

The right hand side of eq. ^ represents work done on the fluid and heating of the 
fluid, respectively. The presence of work and heating terms links disk mechanics with 
thermodynamics. Here, we wish to highlight the dual role of aij-. in eq. 0, it is a term 



- 5 - 



in the transport of mechanical energy flux; coupled in eq. to the strain djVi, it is a 
mechanical energy loss term. In viscous disk models, the latter is, of course, the origin of 
accretion disk luminosity. Since aij is a symmetric tensor, 

{djVi)aij = {l/2){djVi + diVj)aij = {l/2)r]~'^aijaij > (7) 

so the dissipated energy ultimately radiated is necessarily positive definite for incompressible 
turbulence (Landau & Lifschitz 1959). 

In cylindrical coordinates {R, 0, z) the azimuthal equation of motion for an axisymmetric 
viscous disk expresses angular momentum conservation, and takes the explicit form 

dtipRv^) + V ■ (pRv^^^ - riR^Vn) = 0. (8) 
The rotational velocity is Keplerian, 

where M is the central mass, and 

Rn = v^. (10) 

In the simplest form of viscous disk theory, we ignore the vertical structure, treating the 
disk as fiat, set = RQ, and assume axisymmetry. That is, we use the height integrated 
form of the mass, angular momentum, and energy equations. With S denoting the disk 
column density, mass conservation becomes 

dt R OR ' ^ ' 

while angular momentum conservation follows immediately from eq. 

where z/ is the kinematic viscosity, rj = pu. The energy dissipated per unit area of the disk 
Qe ( "emissivity" ) , is found from eq. to be 

Qe = {9/8)um'^. (13) 

Under steady state conditions, the mass flux M = —2TxR'Lvji and the angular 
momentum flux must both be constant. Assuming that the viscous stress vanishes at the 
inner edge of the disk Rq leads to the relation 

m(i-(^)'^') =37rz.S. (14) 

The turbulent paramater z/S is severely restricted by this relation, and can be eliminated 
in favor of M in the expression for the emissivity, leading to (Pringle 1981): 

3GMM / fRQV^\ , , 
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This is the classical (Qe, M) relationship, which leads to a surface emission temperature 
profile Teff{R) oc R~^^^. Its utility lies in the absence of an explicit viscosity term, which 
has been eliminated by the requirements of constant angular momentum and mass flux. 

We may drop the assumption of time steady conditions, forgoing a functional restriction 
for z/S in the process. Using eq. |Tl[] in eq. |]T2[, leads to a generalized mass accretion 
formula: 



T,vrR = 

This in turn may be used back in eq. 
(Lynden-Bell & Pringle 1974): 



oR 



(16) 



1^], yielding an equation for S in terms of u 



'dt 



3_ 9_ 
RdR 



d 



(17) 



This is the classical evolutionary equation commonly used in accretion disk modeling. It 
requires an a priori specification of the functional dependence of v to be useful, and leads 
to diffusive behavior (disk spreading) in eruptive systems. Ultimately mass is transported 
inwards and angular momentum outwards: all of the latter is in a vanishingly small 
component of the former. 



3. Magnetohydrodynamical Turbulence. 



The fundamental assumption underlying essentially all phenomenological modeling 
of turbulent disks is the following: it makes sense to use a two scale approach to 
mathematically represent disk attributes of astrophysical interest. With the possible 
exception of "flickering" in CV systems (e.g., Welsh, Wood, & Horne 1996), observational 
data is assumed to involve length and time scales which are larger than the characteristic 
"eddy turnover" scales of the turbulence. One works with averages that are assumed to 
be well-defined, and to represent the large-scale properties of the disk, much as classical 
dynamo theory (Krause & Radler 1980) represents mean fields and mean helicity. While 
this statement probably does not strike the reader as startling or controversial, it masks 
subtleties. The assumuption has remarkably restrictive consequences. 

The problem is that power spectrum of almost all nondissipative quantities (including 
the stress tensor itself) is dominated by the largest scales of the turbulence, in this case, 
the disk scale height and rotation period. The implicit averaging must be on scales large 
compared to the scale height but small compared with the radius, and large compared with 
the orbital period but small compared with viscous and thermal time scales. There must be 
an asymptotic domain where these scales are cleanly separated, so that a computed radial 
disk profile is insensitive to the averaging procedure. Let us assume this is the case and 
pursue its consequences. 



The dynamical equation of motion in the presence of a magnetic field B is 
p— + {pvV)v = -V \P + — \ -pV$+ I — -V j B + r/yV'v. 
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We have denoted the viscosity as r^y to distinguish it from the resistivity associated with 
the magnetic field, which we will denote as tib- We have dropped terms proportional to 
V-v in the viscous term. The azimuthal component of this equation can be written in a 
form which expresses angular momentum conservation: 



B ( B'^\ 



(19) 



where the subscript p denotes a poloidal vector component. 

We now separate the circular motion RVt from the noncircular motion u, treating the 
latter as a fiuctuating quantity, though not necessarily with vanishing mean. We have 

V = RVte^ + u. (20) 

Although mean drift velocities may be present (the disk must accrete), we assume that such 
motions are small compared with fiuctuations amplitudes, 

(21) 

where the angle brackets denote a suitable average, discussed below. Furthermore, the 
direct contribution to the angular momentum fiux from the microscopic viscosity rjv is 
generally negligible. (This is why turbulence is necessary!) We may drop this term. 

Substituting for v in eq. ([19D , and averaging over azimuth, denoting such means as ()^. 
This gives 

^^{pR^9)^ + V-i? [{pRnup)^ + T] = 0, (22) 

where T is the poloidal stress tensor 

T = {pu^up - B^Bp/An)^ (23) 

We have dropped the pU(f, term in comparison with pRQ in the leading time derivative. 
Henceforth, we will use the Alfven velocity 



ua = (24) 



in favor of the magnetic field vector B. 



3.1. MHD Turbulence and Viscous Disk Theory 



To make contact with the classical viscous disk theory of the previous section we need 
to integrate eq. (^2]) over and to assume that surface terms can be dropped. Furthermore, 
we wish to regardheight-integrated, azimuthal averaged fiow quantities as smooth functions 
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of R. As indicated above, this also implies some sort of radial smoothing — an average over a 
volume small compared with R, but larger than a disk scale height. We follow the notation 
of Balbus & Hawley (1998), and define the density weighted mean of flow attribute X to be 

1 poo fR+AR/2 f2n 

Wp^TT^^/ / . , / pXd(t>dRdz (25) 



27rSAi? j-oo Jr-ar/2 Jo 

where E is the integrated and similarly radially averaged disk column density. The radial 
component of T/p resulting from this operation will be denoted simply as W^i?^. Angular 
momentum conservation becomes 

-(Si?2n), + -— [R'nEiuR), + R'EWr^) = (26) 

Mass conservation follows straightforwardly from integrating the fundamental equation, 
and leads to a form essentially identical to eq. ||11|| : 



9S , 1 d{R^{uR),) 



(27) 



dt R dR 

Using this in eq. (^) gives a general formula for the mass accretion rate 

where the prime ' denotes differentiation with respect to R. Combining eqs. (^) and (^) 
gives us the analogue to eq. (|T7p in a turbulent disk, for any angular momentum profile: 



Since W^i?^ of not known a priori, in practical terms, eq. ( [29| ) represents only a marginal 
improvement on the phenomenological equation (O). But we may see that "viscous" 
evolution does not require the explicit adoption ofa viscous stress tensor. Any disk in 
which ur and (and uar and ma^) are positively correlated must behave similarly, with 
the caveat that the correlation tensor must be a locally defined quantity. 

We have discussed thus far only the dynamics of the turbulence. Once Wr^ is known, 
and it depends primarily on correlations on the largest turbulent scales, the disk evolution 
may be directly calculated by equation (^91) . Classical viscous disk theory also addresses 
the energetics. Since viscosity is the agent of transport, there must be dissipation as well. 
The energy is directly thermalized from its free source in the differential rotation, down to 
thermal scales. 

In a turbulent disk, matters are more complex. Energy cascades from the differential 
rotation to the scales of the largest fluctuations, thence to the integral self-similar scales, 
and finally to the dissipative Kolmogorov scale (which may be set by resistivity rather than 
viscosity). In a steady state disk, we expect that the rate at which energy is extracted 
from the differential rotation, which may be easily calculated in terms of the stress tensor. 
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to be equal to the rate at which it is thermahzed, which would otherwise not be directly 
calculable. The upshot of this is that the steady state turbulent disks behaves viscously in 
their energetics as well in their dynamics. But classical viscous theory makes a stronger 
assumption by its very nature: the rate of thermalization of the free energy of differential 
rotation is the same in both steady and evolutionary models. This may be true in an 
evolving turbulent disk, but it is not obviously true. It depends upon whether the cascade 
is efficient. Fortunately, we shall see that this question may be directly and quantitatively 
answered within the stress tensor formalism we are using. 

The evolution of the magnetic field in a plasma with resistivity is given by 

dB 

— = Vx {vxB -t]bVxB) . (30) 

The energy of mechanical energy equation is obtained by dotting eq. (|I8[) with v, dotting 
eq. ( |30D with B, and combining the two. After some simplification (e.g., Balbus & Hawley 
199877we arrive at 

= PV-v- r]v{d,v,){d,v,) - l^lVxSp (31) 

We have dropped the term — ?7y|V-t;p on the right hand side since in a turbulent disk it 
is generally small compared with the other viscous term. The right hand pressure term, 
though also proportional to V-t;, cannot be dropped; as noted in §2, it represents a link to 
the internal energy of the disk via the first law of thermodynamics. The unwritten energy 
flux in square brackets is 

(1 \ B 

-pv^ + p<^ + Pj +—x{vxB) (32) 

where, as before, we have not included the transport due to explicit viscosity or resistivity. 

The energy flux is more complex than its angular momentum counterpart, but can be 
greatly simplified by retaining only leading terms in a m -C R^l expansion. When averaged 
as before, the radial energy flux is 

s(li?2^2 ^ ^ j:mwR^, (33) 

which may be compared with the radial angular momentum flux of eq. (p6|), 

j:R^n{un)p + ^RWR^. (34) 

(Note that to effect the height integration, we have assumed that the magnetic field is 
force-free above the disk. This is a physically reasonable assumption, but one that is less 
than general. The energy dissipation rate is, however, indifferent to the presence or absence 
of surface terms resulting from vertical integrations, since the ignored vertical fluxes would 
contribute nothing to the local mechanical energy losses in the disk.) The key point is to 
observe that the only turbulence parameters entering into either the angular momentum or 
energy radial fluxes are (u/j)p and Wr^. This is the essence of an a disk. 
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The issue is most clear for a steady model (Balbus & Hawley 1998). In this case the 
accretion rate 

M = -2ttJ:R{ur)p (35) 

is constant, and if the angular momentum flux at the inner edge of the disk Rq is vanishingly 
small, then 



1 



R 



(36) 



Now, one cannot calculate directly the thermalization losses, since the small scale gradients 
are not known. But one can calculate the divergence of the large scale flux, since the spatial 
dependence of Wr^ is determined by angular momentum conservation. This must be the 
small scale dissipation rate. We find 

which is precisely the analogue of eq. ([13|) if IVr^ is replaced by a large scale viscous stress. 
(Note that Qe > 0.) This result can also be obtained directly from the energy equation for 
the u fluctuations themselves, by demanding that sources and sinks balance in steady state 
(Balbus & Hawley 1998). 

More is required, however. In viscous models, the thermalization rate is given by 
eq. (0) whether steady conditions prevail or not. The question before us is whether the 

thermalization rate (ET]) is just as general: does it hold when M is not constant and when 
the energy density oiThe disk changes with time? We now show that it does. 

First, let us recall the fundamental relations for the angular and epicyclic frequencies. 



as well as the specific energy. 



-R^n^ + = —^(R^^). (39) 
2 2RdR^ ' ^ ' 



Thus, the energy flux of eq. (|33D becomes 



^S(««),^(i?2$) + i:RnWn^. (40) 
We may substitute for (u/j)p using equation (^8]). This gives an energy flux of 

" ^C^^^^^'^'^'^ ^ (41) 

The quantity of interest is the heating rate Qe due to turbulent dissipation per unit 
area, and it is given by vertically integrating the left-hand side of eq. (pT]). Making use of 
the above it may be written in the form 



1 dR^^sdT. 1 d 
~ 2R~dR~~dt + rM^ 



(42) 
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If we now use eq. 



Qe 



for 9S / dt in eq. 
1 d 



i?2 



dR 
]_d_ 
RdR 



dR 



we obtain 
R'^Q d 



{R^^yOR 



[^R'Wr^) 



RQ 



{R^y d 



(j^r'^Wr^) - j:R''nwR^ 



(i?3$')' dR 

This unwieldy formula immediately simplifies to 

Furthermore, for any function $, the following identity is easily verified, 

{R-\R^^yy = R-\R^^'y, 

leading to a complete collapse of our expression down to the single term 



R4>) 



Qe = ^W, 



d\nR' 



(43) 



(44) 



(45) 



(46) 



which is the desired result. To leading order in the turbulent fluctuation amplitudes, the 
thermalization rate per unit area of a magnetized disk is given by the above, whether the 
disk is evolving or in a steady state. This result is nicely compatible with classical viscous 
thin disk theory. That not all disk turbulence is so easily subsumed will be seen in the next 
section. 



4. Self-Gravity 



Self gravitational forces can be important for galactic and protostellar disks. In its 
most extreme manifestation, self-gravity can hold the disk together and cause substantial 
deviations from a Keplerian rotation law. But this requires a disk mass comparable to or in 
excess of the central compact mass, and we will not consider this limit. Instead we focus 
on a more common situation in which the local self-gravitating free fall time ~ (Gp)"^/^ is 
comparable to or smaller than the There are several equivalent ways of expressing 

this condition, the classical Toomre (1964) Q criterion being the best known (Binney & 
Tremaine 1987). With cs denoting the sound speed, if 

then local density perturbations are unstable to gravitational collapse in a thin disk. If we 
define the vertical scale height H by cH = Q and the disk density by pH = S, then for a 
Keplerian disk the Q criterion becomes 



(48) 
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in rough agreement with our initial estimate, and we may avoid an exphcit reference to the 
disk temperature. 

Self-gravity is obviously important in the formation stages of a galaxy or a star, but 
it is also likely to be a key component in later evolutionary stages, especially in the outer 
regions of the disk where Vt^ / p is likely to be small. We shall concentrate here on the latter 
case, assuming a well-defined Keplerian disk is present, with self-gravity causing small but 
critical departures from circular flow. The ratio of disk mass to central mass M is found 
from eq. ( ^8|) to be of order (H/QR), so our appoximation is justified for thin disks. 

Progress in the numerical modeling of disk systems has been impressive, and 
sophisticated simulations are now possible, although investigators understandably tend to 
want to explore the more dramatic behavior of very massive disks. One of the interesting 
question these modelers are addressing is whether turbulence wrought by self-gravity is 
amenable to a viscous diffusion treatment (e.g., Laughlin & Royczyska 1996). We now 
examine this point. 



4.1. Dynamical and Energy Fluxes 



The self-gravity potential satisfies the Poisson equation, most conveniently written 
in the form 

P = 4^^^^^*^' (49) 

where the subscript i (or j, k below) denotes a Cartesian coordinate, and the summation 
convention on repeated subscripts is used unless otherwise stated. The connection between 
the gravitational force and its associated stress tensor was first made by Lynden-Bell & 
Kalnajs (1972): 



(50) 



To keep both the gravitational and nongravitational components of the stress tensor on an 
equal footing, define the velocities 



Uc 



(51) 



Then, in the presence of self-gravity and magnetic fields, the R(j) component of the stress 
tensor becomes 

Wr^ = {uru^ + ugrUg^ - uarUa<p)p (52) 

Equations (26) and (29), angular momentum conservation and the disk evolution equation, 
continue to hold in precisely the same form when self-gravity is present, if the stress tensor 
Wr(I, is amended simply as above. Gravitational torques are calculated formally in exactly 
the same way as turbulent and magnetic torques. When Q is of order unity, the kinetic 
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u terms and gravitational uq terms of Wfi^ are comparable if m ~ |V$5|/2f2 — i.e., if the 
fluctuation velocities are due to self-gravity impulses on a rotation time scale. 

Consider next the energetics of self-gravity. Is the volumetric dissipation rate still given 
by equation (46), with the gravitationally amended version of Wr^7 The answer is not in 
general, but under some interesting conditions it is. Let us see how this emerges. 

We seek to write the expression pVidi^s iii conservation form: the time derivative of 
an energy density plus the divergence of a flux. We have. 



(53) 



where the last equality follows from mass conservation plus an integration by parts. There 
is no sign yet of the gravitational stress tensor putting in an appearance, but the final term 
remains dangling for the moment. This may be written 



Returning to vector invariant notation. 



dtlidj^ sM^s)]- (54) 



pvV^c 



9 , ^ 1 



AttG dt 



(55) 



There is some ambiguity as to whether one assigns terms to the energy density or the flux. 
For example, an equivalent formulation of eq. ([55|) is 



pv-V^c 



1 d 



SnCdt 



V$5 (9$5 

AttG dt 



V$i 

SttG dt ^ 



(56) 



But there is no apportionment that of itself produces an R(f) component of the gravitational 
stress tensor. Since the energy flux is most readily interpretable in eq. ([55|), we shall use 
this form of energy conservation in our discussion below. 

The combination — V$5 ■ d^s/9t will be familiar to students of acoustical theory (e.g., 
Lighthill 1978) where precisely this form of energy flux is associated not with a gravitational 
potential, but with the velocity potential of irrotational sound waves. In the acoustic case, 
this emerges from the "Pf " term in the energy flux, a term which is third order in the 
fluctuation amplitudes for incompressible turbulence, and therefore negligible. Indeed, an 
important physical distinction between a disk in which there is a superposition of waves 
and a disk which is truly turbulent is the dominance of the Wpitf, term over the Pv in the 
latter's energy flux. 

If waves were present in a turbulent disk, would this change the relative dominance? 
Not in a thin non-self-gravitating Keplerian disk with good R(j) correlations in the stress 
tensor. In a density wave, the pressure contribution to the energy flux will be of order u'^cs 



in the velocities, whereas the stress tensor term is of order aRQc^. Since u"^ 



ac 



5; 



the 



stress tensor contribution will always be dominant (by a factor of R/H) in a thin disk. 

The appearance of a second order contribution of the potential in the energy flux 
suggests qualitatively new transport features in self-gravitating disks. In retrospect, the 
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breakdown of the a formahsm is perhaps not surprising. Turbulence in hydrodynamical 
shear flows or MHD disks arises because vorticity fields and magnetic fields are "ensnared" 
by shear, and funnel this free energy into fluctuations. These fields may become ensnared 
because both are frozen into their respective fluids. Their evolution is entirely local, and the 
vorticity and magnetic fields are governed by essentially identical equations. Gravitational 
fields are not frozen into the fluid, and we should not expect local dissipation of its 
associated turbulence, which is the inevitable consequence of an energy flux depending 
upon the stress tensor and drift velocity, as may be seen in eq. (33). Self-gravity is generally 
a global phenomenon (its fleld equation is elliptic), and one has no cause to expect a 
repetition of our earlier magnetic success with local theory. 



4.2. The Local Limit 



If instead of the combination pv^s ~ V$5 9t$s appearing in the energy flux, the 
combination Q'V^s ■ d^^s emerged, we would be able to construct a local model of the 
dissipation. In this case, the gravitational component of the stress tensor would couple 
energetically precisely as magnetic and Reynolds stresses couple. It will turn out that the 
vanishing of both these terms corresponds to 




0, 



(57) 



when the disturbances are analyzed in terms of WKB waves. This is a very revealing 
requirement, for it is just this condition that deflnes the corotation resonance in linear 
density wave theory, and it is only at this location that waves couple directly the disk (e.g., 
Goldreich & Tremaine 1979, hereafter GT). It is quite natural, therefore, that this condition 
reemerges as the requirement for gravitationally driven energy stresses to be thermalized. 

Let us examine the structure of the energy conservation equation further. We focus 
for simplicity upon an unmagnetized discs, and assume that RQ S> uq- Denoting the 
volumetric mechanical energy losses as —e, the self-gravitational analogue to equation (1311 
becomes 



d_ 

dt 



$5) + 



8nG 



pi;(y + $ + $5) 



AttG \ dt 



-e. (58) 



We have neglected the pressure contribution to the energy flux (but see below). 



We rewrite the rate of production of mechanical work as given by the left hand side of 
), separating the terms in a suggestive manner: 



d ( 



s 



where the energy density £ is 

£ ■■ 



AnG Dt 



SttG 



(59) 



(60) 
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the local flux 7- is 



and 



+ + (61) 



The final terms on the left side of the eq. (E^) are "anomalous" from the point-of-view of 
a. disk theory, and the fiux will henceforth be denoted as F'^. The radial component of 
is just given by eq. (pS]), with Wji^ amended as in eq. (p^. Were i^^ negligible, we would 
be lead directly to an a disk model via the route we followed in §§ 2 and 3. However, 
since these terms may be of order RQuq in the velocities, they cannot be neglected. Their 
physical interpretation is discussed in the next section. 



4.3. Wave Fluxes in Self-Gravitating Disks 



To understand the role of the anomalous fiux, it is helpful to study it in the 
context of the simplest possible fiuctuating self-gravitating disk model: WKB waves in 
a thin, pressureless disk. (The inclusion of pressure terms leads to a more complicated 
calculation, but with precisely the same final conclusion.) The waves have the canonical 

form exp(z kdx + im0 — iuut), where k is the local radial wavenumber m the azimuthal 
wavenumber variable, and lo the fixed wave frequency, and satisfy the dispersion relation 

(u - mQf = - 2nGJ:\k\. (63) 

The potential $5- has the vertical spatial dependence e~''^'^' out of the disk midplane (Lin & 
Shu 1966, Binney & Tremaine 1987). 

The radial anomalous energy fiux, averaged over azimuth and integrated over height is 

poo 

{P^R^s),dz--—[n-^] I (^^) rf^ (64) 




To do the first integral, we need to be able to express Ur in terms of $5. This relation may 
be read off directly from eq. (11) of GT: 

rriQ — a; , , , , , , , 

= -^i^^siO) sgn(/c). (65) 
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where $5(0) is the midplane {z — 0) value of the potential. Denoting the potential 
amplitude by $5, and assuming p{z) — T:5{z), the first integral is then 

/oo 1 ~ 

^ (puR^s)^ dz = 4^(^^ - sgn(A;) (66) 

The second integral may be evaluated by noting that the integrand depends on z as 
exp(— 2|A;2;|). This gives 



Thus, 



1 



^E'^n- -^("^ - sgn{k) (68) 



The angular momentum flux is also to be found in GT (eq. (30); note their deflnition 
differs by a factor of 27ri? from ours) . It is simply 

^; = -^^|sgn(A;). (69) 

In other words, the anomalous radial energy flux is the product of this angular momentum 
flux and the Doppler-shifted wave pattern speed uj/m — Q. It therefore is identiflablc as 
a true wave energy flux. (Turbulent energy and angular momentum fluxes, by way of 
contrast, are related by a factor Q.) Its significance is that in a Q ~ 1 disk it will contribute 
to the total energy flux at a level comparable to the stress tensor W^i?^ if \uj/m — is 
of order unity. The effect is to prevent self-gravitating disks from behaving like a disks; 
only if the anomalous flux vanishes can a self-gravitating disk behave like a local a disk. 
The "forbidden zone" of wave propagation near the corotation point uj = mQ will display 
properties similar to an a disk. However, when a disk location undergoes forcing due to 
a potential from a wave pattern rotating with a frequency very different from the local 
rotation frequency, it will not behave like an a disk. Such a situation may occur, for 
example, when an exterior disk is forced by the potential caused by a developing central 
bar instability. 

Energy can be exchanged between fluctuations and the differential rotation of the 
disk; unlike angular momentum, it need not be conserved in the noncircular motions. In 
contrast to a turbulent a disk, a self-gravitating disk can evolve by extracting energy from 
the background shear and allocating it to the flucutations (wave energy) without the need 
for mechanical energy dissipation. This allows for angular momentum transport with no 
associated local energy losses. Signiflcant angular momentum transport of this type can 
occur if a global nonaxisymmetric mode develops in an initially gravitationally unstable 
disk. Such a construction need not be merely a transient initial condition. Such features 
are semi-permanent, slowly evolving as the disk background parameters change (Papaloizou 
& Savonije 1991, Papaloizou 1996, Laughlin & Royczyska 1996). More recently, careful 
analyses of self-gravitating disk simulations carried out by Laughlin, Korchagin & Adams 
(1997, 1998) clearly show angular momentum transport produced by the onset of global 
nonaxisymmetric instability and subsequent generation of a global wave pattern which has 
extracted energy from the background shear. Transport of this type, which is seen to persist 
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even after initial saturation, certainly does not have the character of that exhibited by an 
a disk (Laughlin & Royczyska 1996). These simulations, however, involve massive disks 
(comparable stellar and disk masses); similar studies of lower mass Keplerian disks have yet 
to be done. 



4.4. Conditions Under Which Self-Gravity Leads to an a. Disk 



There is a local limit in which the nonlocal energy flux terms vanish and eq. (|37|) is 
recovered. It occurs when the shearing box limit is used to study self-gravitating disks. 

The shearing box approximation is a standard approach to the dynamics of thin disks, 
both self-gravitating (Goldreich & Lynden-Bell 1965, Julian & Toomre 1966, Toomre 1981) 
and non self-gravitating (Goldreich, Goodman, & Narayan 1986; Hawley, Gammie, & 
Balbus 1995). In this scheme, the disk is divided into local Cartesian patches with periodic 
boundary conditions being applied on their boundaries. Thus one considers a small box in 
the disk, and sets up local corotating coordinates corotating with the patch center. Strictly 
periodic boundary conditions are applied in the azimuthal direction, so that 0-averaging 
amounts to averaging over one azimuthal width of the box. However, in the radial direction, 
because of the presence of large scale shear, periodicity is applied at boundary points which 
are azimuthally separating from one another. Related periodic points must shear apart with 
time. Thus, strict periodicity would hold only in comoving, shearing coordinates. Note 
that there is no preferred location in this description so the box may be centered anywhere 
(except, of course, the origin). 

The significance of these boundary conditions is that they force the divergence of i^^ to 
be zero when averaged over the box. If Ai? is the radial extent of the box, and R ± IS.Rj'l 
represent the outer (-I-) and inner (— ) boundaries, then the integrated box average leads to 
a term of the form 



which must vanish. The square bracket notation denotes a difference to be taken between 
the upper and lower indicated locations. The boundary conditions force every fluid element 
on the inner edge to have a corresponding partner on the outer edge, and the appearance 
of D/ Dt, rather than d/dt, ensures cancellation. Were the partial time derivative used, we 
would not be forced to this conclusion, because the disk passes by "faster" at one of the 
boundaries compared with the other. If vanishes (in this averaged sense), we are lead to 
a standard a disk model. Clearly, however, this conclusion is entirely driven by the choice 
of boundary conditions. Energy loss or gain from an evolving wave-like flux would be quite 
incompatible with this type of periodicity. 

There is no physical reason for the above boundary conditions to be satisfied in 
the neighborhood of an arbitrarily chosen disk location. Nevertheless, it is possible that 
there are circumstances under which may in effect vanish. Disks evolving under the 
influence of their own self-gravity tend to hover near the critical Q = 1 level (Laughlin 




(70) 
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& Bodenheimer 1994). WKB waves with radial wavenumber kcrit = vrGS/c| ~ 1/{QH) 
are neutrally stable {ut — mQ = 0); all other wavennumbers propagate. While ostensibly 
comparable in magnitude to may be smaller in a Q = 1 disk. The most responsive 

(dominant?) local modes have u — mQ ^ Q, and this may be enough to suppress F^. 
The effectivenss of this process depends both upon the disk's alacrity in maintaining Q 
near unity, and upon the shape of the wave power spectra. Clearly, numerical simulations 
are needed to resolve the question of whether Q = 1 disks can be treated within the a 
formalism. 

Values of Q near unity are favored, of course, because dropping below this critical 
level results in vigorous dissipative shock heating, raising the temperature and stabilizing. 
Rising above Q = 1 allows the disk to cool and become destabilized (e.g., Sellwood & 
Carlberg 1984). The critical criterion also has some observational support through the 
work of Kennicutt (1989), who has found that the gaseous Q value of active star-forming 
regions of disk galaxies is near critical. More generally, it is yet well understood under 
what conditions heating and cooling will be able to regulate Q efficiently in disks, thereby 
allowing the use of a simplifying a formalism. 



5. Summary 

The dynamical foundations of viscous a disk models are rooted in the correlated 
fluctuations which create the underlying turbulent stresses. In this paper, we have shown 
that the mean flow dynamics of MHD turbulence follows the a prescription, and in 
particular that the disk energy dissipation rate is always give by eq. (^), even if the disk 
is evolving. The local character of MHD disturbances is itself rooted m the flux freezing 
equation, which forces local dissipation of the magnetic field in turbulent flow, analogous to 
vorticity dynamics in an unmagnetized shear layer. 

The mean flow dynamics of a self-gravitating disk in general cannot be described 
so simply. Classical viscous disk theory requires a simple restrictive form for the mean 
momentum and energy fluxes (eqs. and |3^); neither can depend upon transport 



properties other than (u_r)p and Wr^. The energy flux of self-gravitating disks is not 
reducible to a superposition of these quantities. Instead, what we refer to as anomalous 
flux terms are present. These terms allow self-gravitating disturbances (not necessarily of 
WKB form) to propagate nonlocally in the disk via the perturbed gravitational potential; a 
viscous disk cannot communicate with itself in a similar fashion. The angular momentum 
flux (strictly conserved) in a self-gravitating disk has the same canonical form it has in a 
non self-gravitating disk, proportional to Wnf, the energy flux is fundamentally different. 

In a non self-gravitating thin disk, wave energy transport depends upon terms in the 
flux which, while formally present, are small by order H/R. In a self-gravitating disk, the 
additional (non-pressure) terms that are present in the energy flux couple directly to the 
differential rotation of the disk, as does Wr^. This additional coupling means in effect that 
transport becomes global on rotational time scales. Over similar times without self-gravity, 
the domain of wave influence is restricted to the vertical scale-height H. 

Shearing box simulations of self-gravitating disks employ boundary conditions which 
force local behavior, and inevitably must give rise to an a disk. Because self-gravity 
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is intrinsically nonlocal in its manifestations, analyzing transport phenomena in self- 
gravitating disks within the shearing box formalism may be misleading. On the other 
hand, it is possible that critical Q — 1 disks will be dominated by wavenumbers for 
which \uo/m — Vt\/VL is small, in which case a modeling might be a fair phcnomcnological 
description. To date however, global numerical simulations of massive self-gravitating disks, 
do not seem to lend themselves readily to an a formalism. Whether the same is true for 
self-gravitating disks much less massive then their central stars is not yet known. 

We thank J. Hawley for useful discussions. S.A.B. acknowledges support from NASA 
grants NAG5-3058, and NAG5-7500, and NSF grant AST-9423187. J.C.B.P. acknowledges 
support from PPARC Grant GR/H/ 09454. 
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